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C^ ' 1. Introduction 

Models of diffusing particles with hard core interactions were first considered in the 
mathematical literature [1] and the name exclusion process was first coined by Spitzer [2]. 
In the totally asymmetric simple exclusion process (TASEP) particles jump only to the right 
on a one-dimensional lattice but cannot occupy the same site. Mathematical achievements 
include categorising the stationary measures for the process on Z and many results are 
summarised in the books by Liggett [3],ll]. 

Since the early 1990s the TASEP has been of considerable interest within the physics 
community as a prototypical model of nonequilibrium behaviour where, in the steady state, 
a current of particles is supported. In particular the model has been studied on the ring Z^ 



and also on a lattice of length L with open boundary conditions where particles enter at the 
left boundary and leave at the right boundary. Notable achievements have been the use of 
the Bethe ansatz to determine spectral properties of the transition rate matrix [5|, [6], [71 E] 
and the determination of the stationary state in the open boundary case within a matrix 
product formulation [9]. 

A generalisation of the TASEP is to the case of several species of particle. In the 
two-species exclusion process containing first-class particles and second-class particles |T2] 
both first and second-class particles hop to the right with rate 1. However if the site to the 
right of a first-class particle is occupied by a second-class particle the first and second-class 
particle exchange places with rate 1. Thus a second-class particle behaves as a hole from 
the point of view of a first-class particle but behaves as a particle from the point of view 
of a hole. The introduction of such a second class particle is a useful tool to study the 
microscopic structure of shocks [101 El [12] • Besides, the second-class particle problem arises 
naturally from coupling two TASEPs with different densities of particles [I3]: the excess 
particles in the system with higher density acquire the dynamics of second-class particles. 
The stationary state of a system containing second and first-class particles has been obtained 
using the matrix product formulation by Derrida et al. [12]. Based on this work, Ferrari, 
Pontes and Kohayakawa [Tlj introduced a probabilistic construction of the measure. Angel 
[T5] improved this construction providing a combinatorial description of the stationary state. 
In [T6l [T7] , Perrari and Martin showed that Angel's work could be interpreted as a queueing 
system and they generalized it to the A^ species case, for arbitrary A^. 

In the physics literature, the exclusion process with A^ species of particles was considered 
by Mallick, Mallick and Rajewsky [18] and studied for the case A^ = 3. This model, which 
we refer to as the A^- TASEP, is defined by having site variables Tj which may take values 
0, 1 . . . , A^ where A^ is the number of species. (Note that one could alternatively consider 
the state Tj = (a hole) as a species which would imply a total of A^ -|- 1 species; we choose 
instead to use the more common convention.) The dynamics is defined as follows: each bond 
between neighbouring lattice sites has a bell which rings with rate 1. When the bell at bond 
i,i + 1 rings the site variables at i and i + 1 are exchanged provided Tj+i = 0, Tj > or 
Tj+i > Tj > 1. This is equivalent to the following exchanges occurring with rate 1 

a: ^ K for A^ > a: > 1 (1) 

K J ^ J K ioT N> J > K>1 . (2) 

The construction of Perrari and Martin [161 [13 couples A^ realizations of the TASEP in a 
special way, called the A^-line process. To the A^ configurations in the A^-line system one 
associates a configuration of the A^- TASEP. Purthermore, each dynamical event of the A^- 
line process corresponds precisely to a dynamical event in the A^- TASEP. The steady state 
measure of the A^-line system is just a uniform distribution of particles. This implies that 
one may sample the A^- TASEP configurations with their stationary state probability by a 



two step procedure: (a) uniformly sampling a configuration of the A^-line system of particles 
and (b) finding the associated configuration of the A^-TASEP. 

Our aim in this work is to invert this construction to obtain direct expressions for the 
steady state probabilities which generalise those already obtained for the two species case [I2] 
and the three species case [18]. In doing so we shall see how the matrix product formulation 
generalises into a tensor product. 

The paper is structured as follows. In section 2 we review the known results on the 
stationary measure of the 2-TASEP and show how the matrix product representation [12] 
is related to the queueing representation [17]. In Section 3 we consider the A^-TASEP and 
deduce a procedure for computing the stationary state probabilities. In section 4 we construct 
a matrix product representation of the 3-TASEP stationary measure. In section 5 we show 
how matrix product representations of the A^-TASEP may be obtained recursively and we 
conclude in section 6. 

2. Two Species TASEP 

In this section, we review the known solution of the two species TASEP. We also illustrate 
the equivalence between the matrix product solution of Derrida et al., the construction of 
Angel and the queueing process interpretation of Ferrari and Martin. 
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Figure 1. Graphical representation of the construction for two species. 



2.1. Construction of Angel that generates the Stationary State 

We begin by considering the construction of Angel for the two species TASEP on the ring 
Zl- (Note that we often use a different notation to [I5] in order to avoid a clash with some 
standard notation from the matrix product formulation.) The construction is to consider a 
two-line configuration of particles (see Figure [1]). On line 1 there are Pi particles distributed 
randomly (with at most one particle per site) and on line 2 there are Pi + P2 particles 



distributed randomly. Working from right to left we associate to each particle in line 1, the 
nearest particle, at the same site or to the left, in line 2 that has not been associated to 
another particle. The associated particles in line 2 are then labelled 1 and the remaining P2 
unassociated particles are labelled 2. The empty sites of line 2 are labelled thus each site of 
line 2 is labelled 0, 1 or 2. In this way a configuration of the two species TASEP containing 
Pi first-class and P2 second-class particles has been generated through the construction. 
Since we consider periodic boundary conditions, the site at which we begin this procedure 
(chosen as the furthest site to the right in Figured]) does not affect the two species TASEP 
configuration that is generated, but the particular particle in line 2 associated to a given 
particle in line 1 may depend on the initial starting point; for instance, if in Figure[T]one starts 
with the third particle from the right in line 1, then it would be associated to the fourth 
particle from the right in line 2, while the second particle from the right in line 1 would 
be associated to the fifth particle from the right in line 2. As noted in the introduction, 
uniformly sampling the 2-line configurations generates 2- TASEP configurations according to 
their stationary measure. 

Angel showed that by uniformly sampling the two-line configurations, the configurations 
of the two-species TASEP, which we denote C, are sampled with the following probabilities. 

Here oo{Bj) is the weight of the binary string Bj of 0,1 separating the second-class particles 
indexed by j and j + 1 (here j G Zpj indexes the second class particles) . The normalization 

just counts the number of possible 2-line configurations. Note that the form of the measure 
([3]) implies a factorization of the stationary state about the positions of the second-class 
particles. The reason for the factorisation is, as can be seen from Figure [H that all 2- 
line configurations associated with a given 2-TASEP configuration must have the following 
properties: consider a site z, such that there is a particle labelled 2 at i in line 2, then i must 
be empty on line 1; moreover, no particle in line 1 to the right of i can be associated to a 
particle in line 2 to the left of i. This factorization property, which appeared in the matrix 
product formulation of Derrida et al [12], was used in the construction of the stationary 
weights by [T^ . 

The weights uj{B) are given by the following algorithm which we shall refer to as the 
pushing procedure: given the binary string 5, one enumerates the number of strings which 
can be obtained from it by pushing the Is to the right, in addition to the original string. For 
example from the string 110 one obtains 110, 101, Oil. Thus c<;(110) = 3. Similarly, one can 
obtain from 1010 the strings 1010, 0110, 1001, 0101, 0011. Thus o^(lOlO) = 5. 

The measure given by ([3]) is stationary under the dynamics of the 2 species TASEP 
[T2I [m [151 [T9| . Two key properties of this measure are i) the factorisation of the probabilities 



of the 2-TASEP configurations about the position of the second class particles ii) the weights 
uj{B) are given by the pushing procedure described above. 

2.2. Matrix Product Solution of Derrida et al. 

The matrix product formulation has been used to write down the stationary probabilities of 
various interacting particle models, thus allowing models to be solved through the explicit 
computation of physical quantities of interest such as currents, density profiles, correlation 
functions. It was first used to solve the TASEP on a lattice of length L with open boundary 
conditions [9] . It has been extended to the 2-TASEP on the ring Z^^ [12] , partially asymmetric 
processes and more general react ion- diffusion systems (for a review see [20]). In this matrix 
product formulation properties of the stationary measure manifest themselves in algebraic 
relations amongst the matrices involved. Some of these relations have been classified as 
quadratic algebras [2Tj . 

It is important to note that the measure ([3]), along with the calculation of the weights 
uj{B), is equivalent to that first obtained within the matrix product approach, as we now 
show. We recall that we use the variable Tj = 0,1,2 which implies that site i is empty, 
contains a first-class particle or contains a second-class particle, respectively. Let us denote 
by C = (ti, . . . , tl), a configuration of the system. In the matrix product formulation [12] it 
has been proved that the stationary measure may be written as 

P(C) = Z~^W{C) , (5) 

where the weight of the configuration is given by 

L 



W{C) = Tr 



U^r. 



(6) 



and Tr means the trace of the product of matrices X^. The normalization Z = Z{L, Pi, P2) 
is chosen so that the sum of all the probabilities is equal to 1. The matrices X^-. are given by 

Xo = ^, Xi = D, X2 = A, (7) 

that is: if the site is empty we write a matrix E\ if the site contains a first-class particle we 
write a matrix D; if the site contains a second-class particle we write a matrix A. 
The matrices D, E,A obey the algebraic rules 

DE = D + E (8) 

DA = A (9) 

AE = A. (10) 

The only remaining condition to satisfy is that representations of E,D,A may be found which 
give well-defined values for the traces appearing in ([5]). This may be achieved as follows. Let 



In) and {n\ be the column, respectively row, vector having a 1 in the nth coordinate and 
in the other ones, n = 0, 1, 2, . . . . Let A be the projector matrix 



A=|0)(0| 
then D,E may be chosen to be bidiagonal semi-infinite matrices 
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Writing out the matrices explicitly we have 
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where Bj is as in ([3]), and uj{Bj) is now given by 



UJ 



:b,) = {o\i[x^M 



i=l 



(16) 



;i7) 



where / is the length of the binary string Bj and i labels the entries in that string; X^-. is 
either a matrix D or a matrix E according to whether the entry Tj in the string is 1 or 0. 
The algebraic rules (PITU]) imply immediately that the weight of a string comprising a 



segment of consecutive zeros followed by a segment of consecutive ones is equal to 1. In other 
words, a string where any Os are all to the left and any Is are all to the right has weight 1: 



c<;(0---01---l) = {0\E ■ ■ ■ ED ■ ■ ■ D\0) = (0|0) = 1 



'U 



Using rule ([8]) all binary strings can be reduced to strings of the above type and the weight 
of any string is easily computed. For example, 

cu(lO) = (0|D^|0) = (01^10) + {0\E\0) = 2 (19) 

c<;(110) = {0\DDE\0) = (0|DD|0) + {0\DE\0) = 1 + 2 = 3 

c<j(1010) = {0\DEDE\0) = {0\DDE\0) + {0\EDE\0) = 3 + {0\DE\0) = 5 . 

This reduction procedure gives precisely the same result as the pushing procedure of Angel. 
In fact Lemma 2.3 of [14] proves that when u is defined via the pushing procedure, the 
following relation holds 

u;{B10B') = uo{BlB') + uj{B{}B') (20) 

for arbitrary finite binary sequences 5, B' . But this is the same reduction formula that holds 
for the matrix representation: 

cu(SlOfi') = {id\XDEX'\Q) (21) 

= (OIXDX'IO) + (0|XEX'|0) (22) 

where X, X' are the matrix representation of the binary sequences 5, B\ respectively. This 
shows that the definition of uj by the matrix formulation and that by the pushing procedure 
coincide. 

The weights W of 2-TASEP configurations are computed from the weights u of binary 
strings as follows. Recalling that A^-TASEP configurations are translationally invariant under 
the periodic boundary conditions we have: 

iy(0210) = W^(2100) = cj(lOO) = 3, (23) 

1^(0211021) = M/(1021102) = cj(10)a;(110) = 6, (24) 

and the corresponding probabilities are given by 

P(0210) ^ P(2100, ^ |» ^ ^, (25) 

VI/(1021102) 6 , , 

P(02n021) = P(102n02) = 4(^ = ^. (26) 

where Z is defined in (|1]). 

2.3. Queueing Interpretation of Ferrari and Martin 

Ferrari and Martin used A^-line configurations to generate A-TASEP configurations in terms 
of queueing processes. Here, we recall the interpretation in the case of the 2-species TASEP in 
terms of queueing processes and make the connection with the matrix product representation 
of the stationary measure. 

We recall (see e.g. Figure [2]) that a 2-line configuration generates a 2-TASEP 
configuration. Since, as described above, the stationary state factorises about the positions 

7 



Time 



t ^ 



• # 



length of 
the queue 

queue 

arrivals 
services 



Figure 2. Queueing picture for the two species case. Here time t{i) ^ L — i increases from 
right to left, where i is site number. 



of the second-class particles we only need to consider a binary string 5 of Is and Os between 
two 2s in a 2-TASEP configuration. There are several 2-line configurations which generate 
the string B. Those 2-line configurations must satisfy two conditions: (a) both lines 1 and 
2 must contain the same number of particles, equal to the number of Is in B and (b) line 
2 coincides with B. The possibilities for line 1 are then generated from line 2 by pushing 
particles to the right. For example, the lower line of Figure |2] has three strings of type B 
delimited by the three second class particles: Bq = 0, Bi = 1010100 and B2 = 100. The 
upper-line string 1000011 is one of the strings producing Bi, the string 010 is one of the 
strings producing B2 and the string is the only string producing Bq. 

Given a 2-line configuration one can associate to it the trajectory of the length of a 
queue. Consider the labels of the particles of hne 2 as first or second class particles as given 
by Angel's algorithm (illustrated in Figured]). Time for the queue runs from right to left: 
at each site i it is assigned a time t{i) = L — i. The queue has length zero at the times 
corresponding to the positions of second class particles in line 2 (unused service times), a 
particle in line 1 represents an arrival time and a first class particle in line 2 represents a 
service time. At a given time t{i) the length of the queue (constrained to be non-negative) 
increases by one when a particle is present at site i in line 1 but not in line 2 (a new arrival 
occurs and is not serviced); the length of the queue decreases by one when a particle is 
present in line 2 but not in line 1 (a service occurs with no new arrival). If no particles 
are present in lines 1 and 2 (no service or new arrival occurs) or when particles are present 
in both lines 1 and 2 (a new arrival occurs and is serviced) the queue remains at the same 
length. The weight of a 2-TASEP string is then given by all possible queue trajectories, 
consistent with following constraints i) the queue has length zero at the positions of the 
second class particles ii) the effective service times of the queue are fixed by the positions 
of the first class particles. Since the full 2-line configuration can be retrieved by knowing 
the 2-TASEP configuration and the trajectory of the queue, to enumerate the ancestors of a 



2-TASEP configuration it is enough to enumerate the queue trajectories compatible with it. 
We now illustrate how the product of matrices A, D, E defined in flll|12|13l) 
precisely enumerates the possible trajectories of the queue giving rise to a given 2-TASEP 
configuration. The right hand vector |0) represents an initial queue length of 0. At each 
service time of the queue we have a matrix D and at each non- service time a matrix E. A 
vector \n) represents the length of the queue. If the length of the queue just before a service 
time is n > the action of D on In) is 

D\n) = \n) + \n - 1) . (27) 

The two terms represent the two possibilities at the service time: the first represents the 
service of a new arrival at that time, the second represents a service and no new arrival. If 
n = 0, 

D\0) = |0) (28) 

which implies that a new arrival has to be serviced at this time, otherwise there would be 
an unused service which is forbidden. 

Similarly, if the length of the queue is n > just before a non-service time, the action 
of i? on \n) is 

E\n) = \n) + \n + 1) . (29) 

The first term represents no new arrival at that time, the second term represents a new 
arrival at that time. 

The projector |0)(0| at the left end of the string ensures that only trajectories of the 
queue which finish at length are counted and the queue length is set to for the start of 
the next string. 

Remarks 

(i) An alternative way to determine the queue length n at a given time t{i) = L — i is the 
following. In Figure [2] each particle in line 1 is associated to a particle in line 2 by a 
dashed black line. The length of the queue at t{i) is given by the number of dashed 
black lines intersecting a vertical segment passing through i~ , i.e. just to the left of 
site z; vertical dashed black lines do not affect the queue length. 

(ii) The trajectory of the length of the queue in the queueing process with constraints 
described above is precisely a Motzkin path. This allows one to represent matrix product 
calculations in terms of ensembles of Motzkin paths see e.g. [221 1231 [2^ . 
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3. The A^-Species TASEP 

3.1. Construction for N Species of Particle 

In this section we review how the construction for the 2-species case is extended to the A^- 
species case [17]. For A^ species of particle we consider A^-hne configurations of particles. 
The first line comprises Pi particles distributed randomly (with at most one particle per 
site) . The second line comprises Pi + P2 particles distributed randomly and so on until the 
A^th line which comprises Pi + P2 ■ ■ ■ + P/v particles distributed randomly. Initially, in this 
A^-line configuration, particles are not differentiated into species. In the following we define 
the construction by which a species label is attributed to each of the particles. Once this 
has been done the A^th line is identified with an A^- TASEP configuration. 

We start with line 1 and associate each particle in line 1 to a particle in line 2 as in 
Section 12.11 This is done by beginning with a particle in line 1 and associating it with the 
nearest particle, at the same site or to the left, in line 2. We then take the next particle to 
the left in line 1 and associate it in the same way to the first unassociated particle in line 2. 
This process is continued until each of the particles in line 1 is associated with one particle 
in line 2. These Pi particles in line 2 are then labelled 1 and the remaining P2 unassociated 
particles in line 2 are labelled 2. The resulting labels do not depend on which particle we 
began with in line 1, as commented in Section [2m 

We now proceed to associate the particles in line 2 with those in line 3. First we use 
the same procedure as described above to associate the particles labelled 1 in line 2 each 
to a particle in line 3. The Pi associated particles in line 3 are then labelled 1. We then 
proceed to associate the particles labelled 2 in line 2 to P2 of the unassociated particles in 
line 3 (ignoring the particles already labelled 1 in line 3). These particles in line 3 are then 
labelled 2 and the remaining P3 unassociated particles in line 3 are labelled 3. 

The procedure is then repeated up to line A^ and results in P^ of the particles in line 
A^ having label K where 1 < K < N. The construction is illustrated by an example in 
the three species case in Fig. [31 Starting from the random distributions of particles in the 
A^ lines, one obtains a configuration of the A^ species TASEP with Pk particles of species 
K. The probability of the A^- TASEP configuration so obtained is equal to its stationary 
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probability under the A^-TASEP dynamics. This was proven in [T7] ; in Section 14.41 we shall 
give an alternative proof, based on the matrix formulation, for the 3-TASEP. 

Our aim is now to invert this construction. That is, for a given A^-TASEP configuration, 
we wish to compute the probability that it is generated by the above procedure. This amounts 
to the combinatorial problem of counting the number of ways particles may be distributed 
in the A^-line configuration such that the construction will lead to the desired A^-TASEP 
configuration. 

As has been discussed in Section 12.11 for the 2 species case Angel [15] gave such a 
construction and this is equivalent to the matrix product approach of Derrida et al |T2] (see 
section [2l2!) . We showed this by using the queueing representation of Ferrari and Martin. In 
the following we first provide an algorithm, generalising the pushing procedure of Section [2] 
by which the probabilities can be computed. We then construct explicit matrices which 
compute the A^-TASEP weights by book-keeping the generalized pushing procedure. 

3.2. Ferrari and Martin's Construction and the Reverse Algorithm 

As described above, in the A^-species procedure of Ferrari and Martin a configuration of the 
A^-TASEP is obtained from an A^-line configuration, where each of the A^ lines consists of 
a single species TASEP configuration. The procedure may be viewed in the following way: 
from line 1 (a configuration of the single species TASEP) and line 2 one obtains a uniquely 
defined configuration of the 2- TASEP; from that configuration of the 2-TASEP and line 3 
one constructs a configuration of the 3-TASEP and so on, until one reaches a configuration of 
the A^- TASEP. Therefore, a given configuration of the A^- TASEP arises from a whole set of 
(A^ — 1)-TASEP configurations that we shall call its ancestors; each of these (A^— 1)-TASEP 
configurations arises in turn from a whole set of (A^ — 2)-TASEP ancestors etc... 

The stationary weight of the initial A^- TASEP configuration is then given (but for 
an overall normalization constant) by the sum of the weights of the (A^ — 1)-TASEP 
configurations that lead to it (its ancestors). Applying this procedure recursively, we observe 
that this stationary weight is given by the sum of the weights of its (A^ — 2)-TASEP ancestor 
configurations. Finally, because the single species TASEP has a uniform steady state, the 
weight of any A^- TASEP configuration is nothing but the total number of configurations of 
the single species TASEP from which it derives. Therefore, to calculate the weight of a given 
A^- TASEP configuration we must determine the total number of 1-TASEP configurations 
that are its ancestors. 

In the following we shall give a recursive algorithm to determine all the (A^ — 1)-TASEP 
ancestors of a given A^- TASEP configuration. Iterating this algorithm it is possible to obtain 
the total number of 1-TASEP ancestors of a given A^- TASEP configuration; this number 
corresponds to the stationary weight of the A^- TASEP configuration. 

In order to simplify our discussion we shall first present this algorithm for the case of a 
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3-TASEP configuration, i.e. for a string of particles of classes 1, 2, 3 and holes (denoted by 
0). We start from an initial 3-TASEP configuration. 

(i) Freeze the positions of the 2s and the 3s. Construct all possible configurations obtained 
by pushing the Is through the holes towards the right, until they hit a 2 or a 3 (i.e. a 
1 can cross neither a 2 nor a 3). 

— »■ This procedure leads to many 3-TASEP configurations with various positions of the 
Is. From now on the sites where the Is are located will be passive. 

(ii) Keep the positions of the 3s frozen and start moving the 2s. For each configuration 
obtained above, construct all possible configurations obtained by pushing the 2s through 
the holes towards the right, until they hit a 3. Note that the sites occupied by Is are 
spectators and the 2s hop over them as if they do not exist. 

-^ This procedure leads to many 3-TASEP configurations in which the positions of the 
Is and the 2s are fixed. 

(iii) Replace all the 3s by holes. We thus have obtained the complete set of 2-TASEP 
ancestors of the initial 3-TASEP configuration we started with. 

(iv) The stationary weight of the initial 3-TASEP configuration (up to a global normalization 
constant) is the sum of the weights of its 2-TASEP ancestors. 

Let us illustrate the algorithm with the explicit example of the string 2103 

(i) By pushing Is to the right we obtain the strings 2103, 2013 

(ii) By pushing 2s to the right (through the Is), from 2103 we obtain 2103 and 0123 and 
from 2013 we obtain 2013 and 0213 

(iii) We replace 3s by Os to obtain the strings 2100, 0120, 2010, 0210 

(iv) The weight of string 2103 in the 3-TASEP in terms of 2-TASEP weights is given by 

1^(2103) = V^^(2100) + V^^(0120) + ly(2010) + 1^(0210), 

which we may calculate, for example by using the matrix ansatz for the 2-TASEP (15 HTU1) 
and ([23D, as 1^(2103) = 3 + 1 + 2 + 3 = 9. 

It is easy to generalise to the A^-species case and compute the weight of a A^-TASEP 
configuration in terms of the weights of (A^— 1)-TASEP configurations 

(i) Freeze the positions of the species 2, . . . ,K. Construct all possible configurations 
obtained by pushing the Is through the holes towards the right (a 1 cannot cross any 
species of particle). 

(ii) Now in turn for K = 2, . . . , N — 1 push species K to the right keeping the positions 
of species K + 1, . . . ,N frozen and with species 1, . . . , K — 1 spectators, i.e. For each 
configuration obtained from step 1, construct all possible configurations obtained by 
pushing the 2s through the holes to the right, allowing the 2s to hop over Is; then push 
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3s to the right allowing 3s to hop over 2s and Is, and so on until species K — 1 have 
been pushed to the right, hopping over all other species except K. 

(iii) In all of the iV-TASEP configurations generated in step (ii), replace all the A^s by holes. 
We thus have obtained a whole set of (A— 1)-TASEP configurations: this is the complete 
set of ancestors of the initial A-TASEP configuration we started with. 

(iv) The sum of the stationary weight of all these (A^— 1)-TASEP ancestor configurations 
gives the stationary weight of the A^-TASEP configuration. 

We saw that 'the pushing procedure' of Angel is naturally implemented by the D, E and 
A matrices. The algorithm given above is also based on recursive pushing procedures and 
as we shall show in section 14.21 can be encoded by a matrix ansatz; in this case the matrices 
for the A-TASEP are built by using the matrices for the (A— 1)-TASEP as elements. 

3. 3. Queueing Interpretation of N -Species Construction 

Ferrari and Martin also proposed a queueing interpretation for the multiline construction. 
The A^ lines of the A^-species construction correspond to A^ — 1 queues. The first line 
represents arrival times to the queue 1. The second line represents service times for queue 
1. We continue using the convention that the queue time runs from right to left, so that 
the time t{i) associated to site i is given by t{i) = L — i. As we have seen in section 12.31 
unused service times of queue 1 become second-class particles. Then when the particles of 
line 2 have been labelled either first or second-class, they represent the arrivals for queue 2. 
The arrivals are distinguished into first and second-class and the queue is a priority queue: 
at the service times (given by the particles in line 3) the highest priority waiting customer 
is always serviced first. That is, in queue 2 first-class arrivals are served before second-class 
arrivals. The ouput of these service times then become the arrival times for queue 3 with 
unused service times in queue 2 providing third-class arrivals to queue 3. This construction 
is iterated until the particles in line A^ — 1, labelled 1, . . . A — 1 provide the arrivals for 
queue A^ — 1 and the particles in line A^ provide the service times for queue A — 1. When 
the particles in line A^ are labelled 1 ... A they become the output of queue A^ — 1, which 
corresponds to the A^-TASEP configuration. 

4. The Matrix Product Formulation 

In this section, we show how the recursive construction for the A^-TASEP, described in 
Section [3l2| can be encoded within the matrix product formulation, described in Section [2l2l 

4.I. Definition of the Matrix Product Ansatz and Simple Examples 

The matrix ansatz [9j provides a solution to the stationary master equation of the A- 
TASEP (made explicit later in P6|) ). as follows. First consider non-commuting matrices 
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Xq, Xi, . . . , Xjsf, where X^ is associated to particles of class K (in particular Xq is associated 
to holes, Xi is associated to first-class particles etc.). The ansatz represents the stationary 
probability P(C) of a A^-TASEP configuration C = (ri, . . . , tl) (where Tj is equal to K if site 
i is occupied by a particle of class K) as a statistical weight W{C) divided by a normalization 
Z 

P(C) = ^W{C) (30) 

where the weight is given by the trace of the product of L matrices, as follows 

W{C) = Tt{X,,...X,J. (31) 

Here X^-. is equal to Xk if site i is occupied by a particle of class K [K = 0, 1, ... , A^) in 
configuration C. The normalization factor Z (that depends on L and on all the P^'s where 
Pk represents the total number of particles of class K) ensures that Ylc^i^) ~ -'-• ^^ 
emphasize that the matrix formulation depends on the number of species. For example the 
matrices that represents first-class particles in the 2-TASEP and the 3-TASEP are not the 
same. 

If the system contains only first-class particles and holes, it is well known that the 
stationary measure is uniform. Thus the particles and the holes may both be represented by 
one (a scalar) and the matrix ansatz reduces here to a trivial form. 

For the 2-TASEP holes, first-class and second-class particles are represented respectively 
by the matrices Xq = E, Xi = D and X2 = A (in the notation of P, [12]) which satisfy 
f l8|9llT0|) . It is convenient to introduce matrices e and 6 defined by 

E=l + e, D = l + 6 (32) 

where 1 is the identity matrix. Then, by (IHHOfTO!) . the matrices e, A, d generate the following 
quadratic algebra: 

de =1 
dA = 
Ae = . (33) 

4.2. Matrix Ansatz for 3-TASEP 

We now present the matrix product formulation of the stationary state of the 3-TASEP. 

Xi = l®l®D + 5®e®A + 6®l®E (34) 

X2 = A®1®A + A^6®E (35) 

X3 = A^A^E (36) 

Xo = l®l®E+l®e(^A + e®l(^D. (37) 

Note that the matrices are generally sums of tensor products of three semi-infinite matrices 
used in the matrix product representation of the stationary state of the 2-TASEP. From the 
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usual representation of the matrices A, 6 and e given in [9], [12] we obtain explicit expressions 
for the above matrices: they all have a block structure that is bidiagonal. Defining matrices 
F, G, H and R as 
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the matrix representation of (!34] - [37|) reads 
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(38) 



(39) 



(40) 



(41) 



/ 



All these matrices are triply infinite dimensional because their coefficients are themselves 
infinite dimensional matrices with elements D, A and E (which are also infinite dimensional 
matrices). 

4-3. Matrices as Priority Queue Matrices 

We now explain in the case A^ = 3 how these matrices may be obtained from the A^-species 
queueing interpretation of the A^-line configuration discussed in section 13.31 In this case we 
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have a 3-line configuration that represents two queues in tandem. Queue 1 has one type 
of customer which are considered to be first-class: Line 1 represents the arrival times of 
(first-class) customers in queue 1 and line 2 gives the service times of queue 1. The particles 
of line 2 are labelled 1 or 2 according to whether a service time is used or unused. Once 
labelled, the particles of line 2 become the arrival times to queue 2. Queue 2 is a priority 
queue containing first and second-class customers: any first-class customer is served before 
the second-class customers waiting in the queue. The particles of line 3 are the service times 
for queue 2. They are labelled by which class of customer is served; if a service is unused it 
is labelled 3. 

We now consider the possible trajectories of the queue system. To do this we require 3 
integer counters /, m, n: I is the number of first-class customers waiting in queue 2; m is the 
number of second-class customers waiting in queue 2; n is the number of first-class customers 
waiting in queue 1 (i.e. the length of queue 1). The three counters indicate the state of the 
system at each queue time t{i) = L — i (which runs from right to left). The counters are 
therefore indexed by the times t{i) associated to sites i, but we omit this in our notation. 
Remark: The values of the counters can be obtained directly from figure |3] as follows: for 
each site i, the counter I with index t{i) represents the number of black dashed lines crossing 
a vertical segment passing through i~ between lines 2 and 3; m represents the number of 
red dashed lines crossing the same segment and n represents the number of black dashed 
lines crossing the segment between lines 1 and 2; vertical dashed lines are not counted at 
all. The queue counters do not register (a) second class particles served at their arrival time, 
(b) first class particles served in both queues at their arrival time and (c) unused services in 
the third line. However, a 3-TASEP configuration and the trajectories of the three queues 
uniquely determine the 3-line configuration generating it. This implies that it is enough to 
enumerate the set of queues trajectories compatible with the 3-TASEP configuration we are 
computing the weight of. 



The queue counters l,m,n can be represented by a state vector \lmn) 

\lmn) = \l) ® \m) ® \n) (42) 

where |/) = for / < 0. We show now that the matrix product using Xq, Xi, X2, X3 defined 
in ( I3H - I37I) precisely enumerate the possible trajectories of the state of the tandem queues 
giving rise to a given configuration (ri, . . . , tl). 

We list the possible updates of the counters I, m, n at a given site (or time), according to 
the line 3 label of that site, i.e. the site variable Tj in the X-TASEP configuration. Then from 
each possible update of the queue lengths we deduce the necessary action of the matrices Xi, 
i = 0, 1,2,3, on the state vector limn). Finally we can check from the definitions (I3414371) 
of the actions of D,E,A (I27|28II29I) that Xi\lmn) produces the required update of the queue 
counters. 
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Tj = 3 In this case there is an unused service in queue 2 which imphes I = m = 0. In queue 
1 there may or may have not been an arrival therefore n^norn^n+l. Thus, the 
action of X3 must be 

Xs\lmn) = 5ifiSm,o\0) ® |0) O [\n) + \n + 1)] = A®A®E\lmn) (43) 

which recovers the matrix expression for X3, (!36|) . 

Ti = 2 In this case a second-class service occurs in queue 2 which implies that the number of 
first-class customers / = and there is no first-class arrival in queue 2. If there were a 
second-class arrival in queue 2 so that m — > m, it would imply n = as there would have 
to be an unusued service in queue 1. On the other hand, if there were no second-class 
arrival at queue 2 so that m ^ m — 1 then there might or might not be a first-class 
arrival at queue 1 and n^n or n^n + 1. Thus, the action of X2 must be 

X2\lmn) = (5i,o|0)® [5„,,o|"^)® |0) + |m - 1) ® [\n) + \n + I)]] (44) 

= [A(g)l(g) A + A®S(g) E]\lmn) (45) 

which recovers the matrix expression for X2, ( !35l) . 

Tj = 1 In this case a first-class service occurs in queue 2. If there is also a first-class arrival 
at queue 2 then I -^ I, m —* m and n — > n— 1 or n since there is a first-class service 
and possibly a first-class arrival at queue 1. If there is instead a second-class arrival at 
queue 2 (a second-class service in queue 1) then there must be no first-class customers 
in queue 1 and so l^l — l,m^m + l and n = 0. Finally, if there is no arrival to 
queue 2 then there is no departure from queue 1 and there may or may not be an arrival 
at queue 1. Therefore I ^ I — 1, m ^ m and n ^ n or n + 1. 
Thus, the action of Xi must be 

Xi\lmn) = |/)® \m) ® [\n) + \n - 1)] + |/ - 1) ® |m + 1) ® |0)5„,o 
+ |/- 1) ® \m) (g) [|n) + |n+ 1)] 
= |/) ® \m) ® D\n) + 6\l) ® e\m) ® A\n) + 6\l) ® \m) ® E\n) 
= [10 10 D + SiS)eiS) A + 6010 E]\lmn) 

which recovers the matrix expression for Xi, ( l34l) . 

Tj = In this case there is no service at queue 2. If there is first-class arrival at queue 2 
I ^ I + 1, m —>■ m and n ^ n — 1 oi n since there is a first-class service and possibly 
a first-class arrival at queue 1. If there is instead a second-class arrival at queue 2 (a 
second-class service in queue 1) then there must be no first-class customers in queue 1 
and so / — > /, ?7i — i> m -|- 1 and n = 0. Finally, if there is no arrival at queue 2 then 
there is no departure from queue 1 and there may or may not be an arrival at queue 1. 
Therefore I —>■ I, m ^ m and n —>■ n or n + 1. Thus, the action of Xi must be 

Xo\lmn) = |/ + 1) 0\m)0 [|n) + |n-l)] + |/) \m + 1) |0)5„,o 
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+ |/)® \m) ® [|n) + |n + l)] 
= e|/) ® \m) ® D\n) + |/) ® e|m) ® A\n) + |/) (g) |m) ® i?|n) 
= [e^l^D + l^eOA + l^l^^JI/mn) 

which recovers the matrix expression for Xq, ([3] 



4.4- Algebraic Proof of the Matrix Product Ansatz 

The matrix product ansatz may be proved independently of the queueing representation in 
an algebraic way. We shall use the technique of "hat matrices" to prove the ansatz (see e.g., 
[201 [251 [26] for more details). 

We first recall the stationarity condition to be satisfied. The dynamics of the system can 
be encoded in a Markov matrix Q of size fi x fi where Q is the total number of configurations 
of the system. The coefficient Q{C, C) of this matrix represents the rate of transition from a 
configuration C to a different configuration C; —Q{C, C) is the total rate of exit from a given 
configuration C. (Notice that this is the transpose of the usual generator matrix used in 
probability.) Thus the stationary probabilities must satisfy the stationary master equation 

^Q(C,C')P(C') = . (46) 



Due to the local structure of the rules flT|2|) . Q can be written as a sum of local matrices 
that represent the transitions that take place at a bond {i,i + 1) 

L 

QiA+i are (A^ + 1)^ x (A^ + 1)^ matrices whose off diagonal elements Qi,j+i(Tjrj+i; t[t[_^^) give 
the transition rate from configuration rlr'^^i to TjTj+i at the bond z, i + 1, and whose diagonal 
element Qj_j+i(rirj+i; TiTj+i) gives minus the total transition rate out of configuration TiTj+i. 
Since the only transitions involved in the A^-TASEP are exchanges at a bond, we have 

Qi,i+i(JK, KJ) = -Qi,i+i{KJ, KJ) = 1, if K > 1 and J > K or J = , 
Qi^i+i{K'J',KJ) = 0, otherwise, 

where here K, J are indices that take values from to A^. When the steady state probabilities 
are written in the matrix product form (!30!) the local matrix Qi.i+i acts only on the ith and 
the (i + l)th matrices in the product. The stationarity condition (Il6|) then may be written 

L 

^ Tr(X., . . . Xr,_^^,r.^,Xr,^, ... X.J = (48) 



1=1 



where 



K.,r,+i = ^ Qi,i+l{TiTi+i;T-r-^^)Xr^Xr^^^ . (49) 
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rl,rl+l 



That is, 

Ykj = -XkXj for all K > 1 and J > K or J = , 

YjK = XkXj for all K > 1 and J > K or J = , (50) 

Yjj = for all J . 

The key point to prove the validity of the matrix ansatz is to show that l^i,r,+i is a 
divergence-like term, i.e. there exist matrices Xr such that 

^n,n+i = Xr^Xr^^^ - Xr^X^^^^ . (51) 

Summation over i leads to a global cancellation in (148|) . proving thereby that the stationarity 
condition (146!) is satisfied. Combining (I50ll5ip . we obtain the conditions: 

XkXj = XkXj - XkXj for all K > 1 and J > A' or J = , (52) 

XkXj = XjXk - XjXk for all K > 1 and J > K or J = , (53) 

=XjXj-XjXj for all J. (54) 

For A^ = 2, it turns out to be rather easy to solve the above equations (see e.g. |20j): 
indeed, one finds that (152145^ may be satisfied by choosing X^- to be scalars so that they 
commute with X^-. Then fl5^ is immediately satisfied and fl52ll53p reduce to 3 conditions 

^i-'^o = -^1-^0 — -'^o^i (55) 

X1X2 = X2X1 - X1X2 (56) 

X2X0 = X2X0 — ^^^0-^2 ■ (57) 

Choosing Xi = +1, Xq = -1, X2 = and Xi = D, Xo = E, X2 = A recovers §^M)- 

For A^ = 3 it turns out that choosing X^- to be scalars does not allow (l52] - [53l) to be 
satisfied. Thus, the proof rests upon finding the 4 matrices Xk for K = 0, 1, 2, 3. We now 
write explicit forms for the hat matrices that fulfil the above relations when Xk are given 
by (EIHST]): 

li = (1 - (5) ® 1 ® 1 

X2= - A0 6®1 

X3= -A®A®1 

Xo= - Xo + (e - 1) ® 1 ® 1 . (58) 

It remains to verify that relations (I52|53|54l) are satisfied. Here we check a few relations 
involving Xi and Xi. For J = 1 the rhs of ( l54|) becomes 

XiXi - XiXi = {l-6)®l(S)D + 6{l-6)®e®A + 6{l-6)®l®E 

- {1-6)010 D- (1-6)6 ®e0 A -{1-6)6 010 E = . 



Thus fl54j) is satisfied in the case J = 1. 
Using relations ( l33l) we find 

XiXs = A0I0 DA + A060 DE = A0l0A + A060{D + E) , 
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and 

X1X2 - X1X2 = A(g)l(g)A + A^S0E- {-A (g)5(g) D) 

= ^1^2 
X2I1 - X2X1 = A{l-5)(g)l®A + A{l-6)®60E + 

- {-A ®6®D-A6®6e®A-A6^6^E), 



thus fl52ll53p are satisfied for tlie case K = 1, J = 2. Similarly, all relations fl52|53ll54l) may 
be verified. 

5. Hierarchical Matrix Ansatz for the Multispecies ASEP 

In this section, we generalize the previous construction to the multispecies totally asymmetric 
exclusion process on the ring Z^, with N classes of particles for any A^ > 1. We show that a 
matrix ansatz for a system containing A^ classes of particles (plus holes) can be constructed 
recursively knowing a matrix ansatz for a system with A^ — 1 classes of particles (plus holes). 
We shall simply present the results here and give some examples; we leave the algebraic 
proof and further generalizations to a future publication. 

The matrices Xj^ at level A^ are obtained by making tensor products of the X|,j ~ 
defined at level A^ — 1 with some matrices axAi constructed from e, A, 6 and 1. The matrix 
ansatz is given by 

7V-1 
'W _ V^ (N) „ v(N~l) 



M=0 

N-1 

Xf ) = aifo) ® Xf "^) + J2 «S ® 4f "'^ for 1 < K < A^ . (60) 



M=K 



We emphasize that in this section our notation for the matrix Xj^ has two indices: the 
lower index K denotes the class of the particle represented by the matrix, whereas the upper 
index A^ gives the total number of classes considered in the system. 

The fundamental building blocks to construct the a)^l^ matrices are the matrices e, A, 6 
and 1 (identity). The a^^ ^^^ then given by 



a\;'^> = i®(^-i) 






(61) 


^{Nj ^ ^«(M-i) ^ ^ ^ ^m-M-i) for 1 < Af < AT - 1 . 






(62) 


For X > 1 we have 








^{N) ^ ^^(K-l) ^ g ^ -^m-K-1) ioil<K<N-l 






(63) 


a^^l = A«(^-i) ® 1^(^-^) 






(64) 


afl^ = A®(^-i) ® 6 ® l^(*^-^-i) ® e ® l®(^-^-^-i) for 1 < K 


< 


M <N 


-1(65) 
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AN) _ 4®(7V-1) 
"atQ ~ ^ 

It is understood in the formulae above that any matrix raised to a tensor-power equal to 
zero is equal to the scalar 1 which can be removed from the tensor product. 

Note from fl59|60l) that the matrices X^ at level A^ are composed of tensor products 
of (2) fundamental matrices e,A,6 or 1. 

5.1. Some Examples 

Using the hierarchical matrix ansatz given above, we study explicitly the cases A^ < 3. 

For A^ = 0, the system does not contain any particles but only holes. There is only one 
configuration which has probability 1. Thus, we define Xq = 1. 

For A^ = 1, we obtain from equations ( l59l) and (160|) . using the fact that Xq = 1 

Xi') = aS and xf ^ = a^^ . (67) 

But from equations (16T|) and (!66|) we find that Oqq = a^Q = 1 and we recover the fact that 
for a system with only one class of particles the stationary measure is uniform and therefore 
the matrix ansatz is trivial. 

For X = 2, we find from equations (EI]), ([62]), (El) and ([66)), that 4J = 1, 4? = e, 

A. Then, from the recursion relations ( l59l) and ( l60l) . and using 
3 that 

= E, Xf ^ = ajj + aS? = 1 + 6 = D, (68) 

(69) 



alfj = 6, a[^ = 1, and ago 
the fact that X^^^ = x[^^ -- 


= A T 

= 1, we 


^her 
dec 








Y^ _ A-i) 
^0 ~ "00 

and Xf ^ = 


+ 4? = 

- "20 ~ 


= 1 

A. 



We retrieve the fundamental matrix ansatz ( l32l) . 

For X = 3, using the recursion relations (!59l) and (160!) . the definitions (16T1466I) and the 
results (l68l) ■ fl69l) . we obtain 



(70) 



(71) 

L ® A (72) 

(73) 

and we retrieve the expressions ( l34] - [37l) . 

6. Discussion 

In this work we have considered the multispecies totally asymmetric exclusion process on 
the ring Z^, (although our results are generalisable to Z). We have shown how the stationary 

21 



-(3) 



= ag ® X(^) + 4? ® Xf ) + 4? ® Xf ) 




= l®l®E + e®l®D + l®e®A 


-(3) 
1 


= aS?®xf +af,)®X;naS®xf) 




=5®l®E+l®l®D+5®e®A 


-(3) 
2 


= ag ® Xf^ + ag ®xf^ = A®5®E + A 


-(3) 
3 


= affl ®xl^^^ =A®A®E, 



measure may be written in a matrix product formulation, thus providing an algebraic proof 
of the stationary measure which we presented for the three species case A^ = 3. For 
arbitrary A^ we have shown how the matrix product formulation may be constructed in 
a hierarchical fashion, although we leave the algebraic proof of the stationary measure to a 
future publication. 

Ferrari and Martin have constructed the stationary state of the A^-TASEP as the output 
of A^ queues in series with A^ priority-classes of customers, for all A^. The construction takes 
a A^-line binary configuration sampled at random and produces a A^-TASEP configuration 
whose resulting law is invariant for the A^-TASEP. We have shown that the matrix ansatz for 
A^ = 2 of Derrida et al [12] gives a mechanism to count the number of 2-line configurations 
producing a given 2-TASEP configuration; the matrices may be thought of as acting on the 
space of queue counters. For N = 2 there is just a single queue with one type of customer 
and the queue counter is simply the length of the queue. We have also extended the matrix 
ansatz for A^ > 2. In this case there are multiple priority queues and there are several queue 
counters representing the number of each class of customer in each queue. This results in the 
queue matrices acting on tensor product spaces and accordingly the 'matrices' of the matrix 
formulation become higher rank tensors. This relation between the matrices and queueing 
processes also provides us with a natural representation of the space on which the matrices 
act. Until now, it was believed that the matrices act on a purely formal 'auxiliary' space 
which did not have any physical interpretation. 

The algebraic proof of the stationary measure for A^ > 3 relies on the existence of 'hat' 
matrices [251 1261 [20] described in Section 14. 4[ This is in contrast to the N = 2 case where 
the hat matrices were simply scalars and the relations obeyed by the matrices become a 
quadratic algebra, as in the A^ = 1 open-boundaries case of [9]. For N > 2 the relations 
between the matrices have a more complicated algebraic structure and it would be of interest 
to explore this further. 

One advantage of the matrix product formulation of the stationary measure is that it 
provides a framework within which the calculation of quantities of physical interest, such as 
correlation functions, can be carried out. So far we have not attempted such calculations 
but it would be important to do so. 

Finally, we mention that the multispecies TASEP may be generalised in several ways 
by introducing rates which differ from one or additional processes. For example, for A^ = 1 
allowing particles to carry out forward exchanges with holes with rate p and backward 
exchanges with rate q generates the partially asymmetric exclusion process for which a 
matrix product formulation of the steady state on the open boundary system has been fully 
worked out [271 128] • In the case N = 2 there are partially asymmetric generalisations which 
admit matrix product stationary states [121 EHj • So far, in work in progress, we have found 
a partially asymmetric generalization of the matrix product ansatz presented in Section [51 
It would be of interest to understand how the queueing interpretation of the steady state 
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should be modified. 
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